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Abstract 

We compute the 0{g'^) contribution to the thermodynamic pressure for Wil- 
son fermions in the standard, the twisted mass, and cfover improved formulation 
in lattice perturbation theory, including finite mass effects. We compare the con- 
tinuum approaches of these discretizations for the massive ideal and interacting 
gas. In all cases, for > 8 cutoff effects of Wilson type fermions are comparable 
to those of staggered fermions, but asymptotic scaling requires Nr > 10. 



1 Introduction 

The most prominent systematic errors in numerical studies of lattice QCD are due to 
cutoff effects. The goal of improvement is to remove the leading lattice corrections 
for specified observables in a series expansion in the lattice spacing a, either by lattice 
perturbation theory or by nonperturbative means [H |2]. The twisted mass fermion 
formulation represents a special case, being automatically 0{a) improved at maximal 
twist [3]. Corresponding nonperturbative quenched studies are reported in [1]. 

For studies of the QCD equation of state at finite temperatures, improvement is 
particularly relevant because a large fraction of the thermally distributed field modes 
lives on the scale ~ T = I/qNt-, i.e. on the scale of the cutoff. In perturbation 
theory, the pressure can be expressed as p = p^^ + p^^ + . . . + pp^ + p^P + . . . where 
p'q^ and p^p^ denote the pure gauge and fermionic contributions to the order of the 
coupling (yf", respectively. Despite existing three- loop results for the free energy at zero 
temperature [5], perturbative studies of Wilson-type fermions at finite temperature 
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are restricted to the ideal gas and chiral limits pp (m = 0); cf. [6] and references 
therein. However, these limits are far from the phenomenologically relevant situation, 
T < 4Tc, and with quark masses still heavier than physical in many cases (cf. |7] for a 
recent review). For twisted mass fermions there is an investigation only at vanishing 
temperature and at tree level [8]. 

In this paper, we evaluate the effects of nonzero quark masses and interactions to 
order 0{g'^) on the size of cutoff effects in the thermodynamic pressure by means of 
lattice perturbation theory. It is well known that the leading cutoff effects of inter- 
acting staggered fermions to the pressure are of O(a^), whereas massive, interacting 
Wilson fermions in general will have 0{a) effects. Clover improved Wilson fermions 
and twisted mass fermions both lead to an 0{a) improvement if the respective param- 
eters are tuned to a suitable value. Here, we explicitly evaluate the next-to-leading 
order (NLO) contribution to the pressure for these discretizations, and quantitatively 
compare their cutoff effects with previous results for staggered fermions [9]. We find 
that the interactions increase cutoff effects and thus should not be neglected for making 
an optimal choice of improvement scheme. To leading order in the interaction and for 
lattices with Nr > 8, Wilson fermions are competitive with staggered fermions. 



2 Wilson fermion actions 

The twisted mass formulation is obtained from the standard Wilson action Sw for 
Nf = 2 degenerate fermions as 

5'tm = Swim) +i/x ^ a'^^/jf{x)-f5Tfj,i/jf{x), (1) 
2;,/,/' 

where we will always work with the standard choice of the Wilson parameter, r = 1. 

denotes the twisted mass parameter and the diagonal Pauli matrix acts in flavour 
space. The bare quark mass in this formulation is rriq = \/vp? + /i^, 

m = niq cos(a;), fi = niq sin(ci;). (2) 

Twisted mass fermions provide an automatic 0{a) improvement over Wilson fermions 
at maximal twist, u = 7i/2 (see [IQj). The free quark propagator for twisted mass 
fermions is 

^tm(p) = a -— -2 -— , (3) 

p + [2P + am) + (a/i)^ 

with the dimensionless lattice momenta p^ = sin(ap^) and p^ = 2sin(ap^/2). 
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To identify cutoff effects, the dispersion relation E{p), which is obtained from the 
poles of the propagator, can be expanded in small lattice spacing, 

E{p = 0) = mg- ^aml cos{uj) + 0{a^). (4) 

Cutoff effects of 0{a) set in with quark mass, and their removal at maximal twist is 
apparent. 

An improvement scheme that is frequently used in numerical simulations, at both 
zero and finite temperatures, are clover improved Wilson fermions [TT] . 

Ssw = Sw{m) + Csw iga^ ^ (5) 

with the gluon field strength -F^J^. In this case one uses the propagator of standard 
Wilson fermions, but with rescaled bare quark mass and gauge coupling [121 131 [S], 

mf^ = {m-mc){l + bmam), (6) 
9sw = a^i^ + bg arriq). 

Here mc{g'^) denotes the additive quark mass shift corresponding to the chiral limit, 
which as &m(fl'^) and bg{g^) can be determined in perturbation theory. The tree- level 
dispersion relation is 0{a) improved. 
The coefficient 

Csw can also be determined in perturbation theory. For our purposes, 
we need these various coefficients to O^g'^) (cgw, bg) or 0{g'^) [HI [TSl [IS]- Since these 
quantities are determined to improve the dispersion relations entering the propagator, 
they can be applied to the finite temperature case as well as to zero temperature. 

3 The ideal gas limit 

The thermodynamic pressure is determined by the logarithm of the partition function 
as p = {(3V)~^lnZ, i.e. the sum of all bubble diagrams. The ideal gas contribution 
of Nf = 2 free Wilson or twisted mass fermions is immediately obtained from the 
logarithm of the inverse propagator, where as usual the divergent vacuum pressure has 
to be subtracted by way of renormalization. In order to evaluate the one-loop integrals 
numerically, we compute the corresponding sums for finite volume explicitly, finding 
that a smooth extrapolation to the thermodynamic limit can be done. 

Cutoff effects for the free gas are obtained by comparison with the continuum Stefan- 
Boltzmann or ideal gas limit, which contains a one- dimensional integral that has to be 
solved numerically in case of finite mass (cf. also |17]). 
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Figure 1: Left: Stefan-Boltzmann limit for massless fermions, normalized to the contin- 
uum result. The solid line is the analytic (9(l/iV^) prediction. Right: Mass dependence 
of the Stefan-Boltzmann limit on a lattice with Nt- = 8. 
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Figure 2: Left: Mass dependence of the 0{a^) contribution to the pressure. Right: 
Comparison of the first coefficients for the pressure series in l/Nr- 



In [6] the lattice pressure for quarks has been expanded into a power series in I/Nt-: 

rp4 2 



P_ _ [oN^ J^3p ^ g-iV.a£(p,H) (7) 

p^"^) 1 I 1 



— P _|_ ' 1 _|_ P 1 _|- (S.\ 

rp4 \ rpi • • ■ \°) 

For massless quarks, the coefficients p'^°'^\m) /T^ and p^"'*\m)/T'^ can be calculated in 
closed form [6] and remarkably are the same for standard staggered and Wilson-type 
fermions. Differences are introduced only in higher orders of the lattice spacing. We 
have repeated this procedure for massive fermions. In this case, one is left with one- 
dimensional integrals for the expansion coefficients depending on the quark mass that 
can be solved easily by numerical integration. 

The fermionic contribution to the pressure on the lattice is shown normalized to the 
continuum contribution in Fig. [1] The chiral limit has been discussed before [6] and is 
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reproduced on the left for comparison. In this case there is no difference between the 
Wilson formulations; cf. Eqs. ( 12|6|) . Even though in the chiral limit the leading cutoff 
effects are found to be identical for both types of discretization, staggered fermions 
show a more rapid continuum approach than Wilson fermions. As expected from 
the dispersion relation, Eq. (jl]), switching on a mass increases the cutoff effects for 
Wilson fermions to 0{a), whereas massive staggered fermions keep scaling as (9(a^). 
Fig. [1] (right) shows that cutoff effects are generally increased by finite quark masses. 
However, the effect becomes sizeable only for large masses and can be alleviated by 
0{a) improvement. 

The analytical approach supplements the numerical findings. Up to large quark 
masses the differences between the 0(a.^) contributions for the fermion formulations 
that have been considered are very small (see Fig. [2], left). Furthermore, the 0{a) 
contribution to pure Wilson fermions is almost negligible when compared to the O(a^) 
(see Fig [21 right). This will only change only once > 100. 



The leading 0{g'^) corrections to the fermionic pressure due to interactions are given 
by two diagrams. 



Twisted mass fermions and ordinary Wilson fermions can be treated simultaneously 
using the propagator Eq. ([3]), with identical vertices. We use the expressions as given 
in [2] with the gauge propagator from the standard Wilson plaquette action. 

For purposes of comparison with staggered fermions, we use the corresponding result 
given in the appendix of [9]. All integrals have been evaluated numerically using 
integration routines of the CuBA-library [18]. Finally, in order to extract the cutoff 
effects we also need the continuum 0(5f^)-corrections from [191 120] . 

Unfortunately, the difference between the finite A^^ and the vacuum contributions, 
which is the quantity of interest, shrinks rapidly ~ 1/A^^ and for A^^ = 8 is only about 
6% of the numerically evaluated integrals, rendering an accurate evaluation difficult. 

Fig.|3]shows a comparison of the NLO contribution for Wilson and staggered fermions 
in the chiral limit. As in the noninteracting case, staggered fermions have smaller cutoff 
effects on coarse lattices. Comparing the scales of Fig. [3] and Fig. [T] (left), interactions 
indeed appear to make this difference more pronounced. However, also for staggered 
fermions a? scaling is not setting in until At > 8 at least, by which point the cutoff 
effects of the Wilson discretizations are comparable. 



4 Weak coupling expansion to 0{g^) 
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Figure 3: Two-loop contribution to the pressure for massless fermions normalized to 
the continuum result. 
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Figure 4: Continuum approach of the pressure of massive fermions with a fixed renor- 
malized quark mass vn^jT = 0.03. Left: Ideal gas. Right: Two-loop contribution. 



5 Continuum limit for constant mass 

For finite mass, mass renormalization leads to an additional contribution to 0{g'^) from 
the one-loop integrals. For all quark masses, we have = Z^mq. For Wilson fermions 
there is additive quark mass renormalization = mo—mc, while for both the staggered 
and the twisted masses there is only multiplicative renormalization. The renormaliza- 
tion constants = 1 + Zm g"^ + . . . are in principle discretization dependent. How- 
ever, the needed first nontrivial coefficient is universal [H], Zm = —QCp ln(a/is)/(47r)^. 
Choosing the scale fig ~ 2ttT, we find that the scale dependence is small in all cases. 

The behavior of massive fermions is shown in Fig. HI Despite the rather large numeri- 
cal uncertainties, one finds a qualitative correspondence to the one-loop case, especially 
the small dependence on the quark mass and the particular discretization. As in the 
free case, on lattices A^,- > 8 Wilson fermions are found to be competitive with stag- 
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gered fermions. However, a? scaling is not observed for any of the discretizations before 
^ 10. Up to that, 0{a) effects for unimproved Wilson fermions are negligible. The 
comparison of the scales between the free and the interacting situation shows that con- 
siderations based on the ideal gas limit alone may be insufficient for the discussion of 
improvement schemes. Note that the absolute size of the cutoff effects for the unim- 
proved action appears smaller than for the clover action. This indicates significant 
higher order contributions for the considered here. 

6 Discussion and conclusions 

Our results indicate that the statement made in [6] for the massless ideal gas limit, 
namely, that an improved dispersion relation leads to the same level of improvement in 
the pressure, might generalize to the massive and interacting case. Let us now discuss 
the relevance of our results for numerical simulations. In a regime, where perturbation 
theory is valid, we can consider the relative cutoff error for the pressure, 

I^PfI /PF.cont = \VFM - PF,cont| /PF.cont, (10) 

where pF,cont and pF,\at are each calculated to some specified order in the coupling, in 
our case 0{g^). In the limit of weak couplings the denominator is dominated by the 
ideal gas limit, and the cutoff effects of the interactions get normalized to that number, 
rather than to p^T^conf ^^i^ ^^^^ they play a rather unimportant role quantitatively, 
for all discretizations. In fact, the sign of the leading order corrections is negative, so 
that the relative cutoff error actually begins to shrink when the coupling is increased 
from zero. However, one cannot draw any conclusions from this for the behavior of 
the pressure near Tc. As the coupling grows, Eq. flTU]) quickly develops a pole beyond 
which results cannot be extrapolated. For longer series the pole might disappear, but 
the relative cutoff error as a rational function will be nonmonotonic in general. 

Hence, all one can say in a perturbative analysis is that the lower the temperature, 
the higher the relative importance of Piriat/PFcont- Therefore our analysis suggests 
that for simulations near T^, Wilson fermions scale comparably to standard staggered 
fermions for fine enough lattices, N-r > 8. 

We conclude that interactions play a crucial role in discussing improvement schemes 
for lattice QCD simulations at finite temperatures. An advantage of maximally twisted 
fermions in this context is that its 0{a) improvement appears to hold also nonper- 
turbatively. This suggests to further explore the use of this discretization also for 
finite temperature simulations [21]. However, independent of the improvement scheme 
chosen, a^-scaling appears to set in only on lattices A^,- > 10. 
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